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Abstract 

This paper revisits a well-studied anti-plane shear deformation problem formulated 
by Knowles in 1976 and analytical solutions in general nonlinear elasticity proposed by 
Gao since 1998. Based on minimum potential principle, a well-determined fully nonlinear 
system is obtained for isochoric deformation, which admits non-trivial states of finite anti¬ 
plane shear without ellipticity constraint. By using canonical duality theory, a complete 
set of analytical solutions are obtained for 3-D finite deformation problems governed by 
generalized neo-Hookean model. Both global and local extremal solutions to the noncon- 
vex variational problem are identified by a triality theory. Connection between challenges 
in nonconvex analysis and NP-hard problems in computational science is revealed. It is 
proved that the ellipticity condition for general fully nonlinear boundary value problems 
depends not only on differential operators, but also sensitively on the external force field. 

The homogenous hyper-elasticity for general anti-plane shear deformation must be gov¬ 
erned by the generalized neo-Hookean model. Knowles’ over-determined system is simply 
due to a pseudo-Lagrange multiplier and two extra equilibrium conditions in the plane. 

The constitutive condition in his theorems is naturally satisfied with b = A/2. His ellip¬ 
ticity condition is neither necessary nor sufficient for general homogeneous materials to 
admit nontrivial states of anti-plane shear. 
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1 Remarks on Nonconvex Variational Problem and Challenges 

Minimum total potential energy principle plays a fundamental role in continuum mechan¬ 
ics, especially for hyper-elasticity. One important feature is that the equilibrium equations 
obtained (under certain regularity conditions) by this principle are naturally compatible. 
Therefore, instead of the local method adopted by Knowles [HIE], the discussion of this 
paper begins from the minimum potential variational problem ((V) for short): 

(V) : min jll(x) = jf W(V X )dB - ^ X ■ tdS] X € Y c | , (1) 

where the unknown deformation x(x) = { X i(xj)} S T a is a vector-valued mapping B C R 3 — > 
ui C M 3 from a given material particle x = { Xi } £ B in the undeformed body to a position 
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vector in the deformed configuration ui. The body is fixed on the boundary S x C dB , while 
on the remaining boundary St = S x n dB , the body is subjected to a given surface traction 
t(x). In this paper, we assume 

* 0 = £ >V 1 ’ 1 (i5;M 3 )| det(Vx) > 0 x(x) = 0 Vx e S x }, (2) 

where W 1,1 is the standard notation for Sobolev space, i.e. a function space in which both 
X and its weak derivative V% have a finite T 1 (I3) norm. Clearly, a function in W 1,1 is not 
necessarily to be smooth, or even continuous. For homogeneous hyperelastic body, the strain 
energy W(F) is assumed to be C 1 on its domain F a C M 3x3 , in which certain necessary 
constitutive constraints are included, such as 

W(F) >0 VFG F a , W'(F) ^ocas ||F|| -> oo. (3) 


For incompressible materials, the condition det F = 1 should be included. Finally, X c = 
{x £ X a \ V% £ F a } is the kinetically admissible space, which is nonconvex due to nonlinear 
constraints such as det(Vx) > 0. Also, the stored energy W(F) is in general nonconvex. 
Therefore, the nonconvex variational problem ( V ) has usually multiple local optimal solutions. 
Let Aj, C X c be a subspace with two additional conditions 

X b = { X €X c \ x£C 2 (B;R 3 ), W(F(x)) £ C 2 (F a ;R)}, (4) 


the criticality condition <5II(x; $x) = 0 V5 x £ leads to a nonlinear boundary-value problem 


(BVP) : 


f —V • er(Vx) = 0 in B, 

\ N • <r(V%) = t on St, X = 0 on S x 


(5) 


where, N 6 K 3 is a unit vector normal to OB, and <r(F) is the first Piola-Kirchhoff stress 
(force per unit undeformed area), defined by 


<r = VW(F), or ^ = i,j = l,2,3 


dFi 


( 6 ) 
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which is also a two-point tensor. 


Remark 1 (KKT Conditions, Isochoric Deformation, pseudo-Lagrange Multiplier) 

Strictly speaking, there is an inequality constraint in X c , i.e. the admissible deformation con¬ 
dition det(Vx) > 0. According to the mathematical theory of non-monotone variational 
inequality,in addition to the equilibrium equations in (BVP), we have the following KKT 
conditions 

P < 0, det(Vx) > 0, pdet(Vx) = 0 (7) 

where p is a Lagrange multiplier and p < 0 is called the condition of constraint qualification. 
The equality pdet(V%) = 0 is the well-known complementarity condition in variational in¬ 
equality theory, by which we must have p = 0 in order to guarantee the inequality constraint 
det(V%) > 0. Therefore, this constraint is actually not active to the problem (V). Such an 
inactive constraint is not a variational constraint. 

For incompressible deformation, the inequality condition det(V%) > 0 in X c should be 
replaced by an equality constraint det(V%) = 1. In this case, (V) is a constrained variational 
problem. The KKT conditions ([7]) should be replaced by (see [16] ) 

P + 0, det F(x) = 1, p(det F(x) - 1) = 0. (8) 
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and we must have p(x) / 0 for a.e. x e B in order to ensure det F(x) — 1 = 0. The associated 
{BVP) should be 


(RVP) ■ I ~ V ' = °> det(Vx) = 1 in B, 

1 >p ’ \ N • cr(Vx,p) = t on 5 t , X = 0 on S x . 


(9) 


in which, a{F,p) = VIT(F) — pF~ T , where F~ T = (F T ) _1 . In this case, we have two variables 
(x,p) and two equations in B , thus, the problem ( BVP) p is a well-defined system. 

For isochoric (i.e. volume preserving) deformation, say the anti-plane shear problems, the 
condition det F = 1 is trivially satisfied and the complementarity condition p(det F — 1) = 
0 Vp(x) ^ 0 a.e. in B. In this case, the trivial condition detF = 1 is not a variational 
constraint for ( V ) and the arbitrary parameter p(x) is not an unknown variable for ( BVP) p . 
Otherwise, the ( BVP) p is an over-determined system. This fact in KKT theory is important 
for understanding Knowles’ anti-plane shear problem. Such a parameter for trivial condition 
can be called pseudo-Lagrange multiplier. 

Physically speaking, the hydrostatic pressure p is not necessary to be zero even for isochoric 
deformations. There are many examples in the literature, see the celebrated book by Ogden 
EBj as well as many famous papers by Rivlin on volume-preserving deformations of isotropic 
materials (simple shear, torsion, flexure, etc.)[^] ■ 


Remark 2 (Convexity, Multi-Solutions, and NP-Hard Problems) 

The stored energy W[F) in nonlinear elasticity is generally nonconvex. It turns out that 
the fully nonlinear {BVP) could have multiple solutions {x fc (x)| at each material point 
x G B s C B. As long as the continuous domain B s / 0, this solution set {x fc (x)j can form 
infinitely many solutions to (BVP) even B C M. It is impossible to use traditional con¬ 
vexity and ellipticity conditions to identify global minimizer among all these local solutions. 
Gao and Ogden discovered in |10j that for certain given external force field, both global and 
local extremum solutions are nonsmooth and can’t be obtained by Newton-type numerical 
methods. Therefore, Problem (V) is much more difficult than {BVP). In computational 
mechanics, any direct numerical method for solving {V) will lead to a nonconvex minimiza¬ 
tion problem. Due to the lack of global optimality condition, it is a well-known challenging 
task to solve nonconvex minimization problems by traditional methods. Therefore, in com¬ 
putational sciences most nonconvex minimization problems are considered to be NP-hard 
(Non-deterministic Polynomial-time hard) [llj. 

Direct methods for solving nonconvex variational problems in finite elasticity have been 
studies extensively during the last fifty years and many generalized convexities, such as poly-, 
quasi- and rank-one convexities, have been proposed. For a given function W : J-" a M, the 
following statements are well-known (see [23|)[^J 

convex => poly-convex => quasi-convex => rank-one convex. 


Although the generalized convexities have been well-studied for general function W{ F) 
on matrix space M mxn , these mathematical concepts provide only necessary conditions for 
local minimal solutions, and can’t be applied to general finite deformation problems. In 
reality, the stored energy W{F) must be nonconvex in order to model real-world phenomena, 

1 Personal communications with David Steigmann, Ray Ogden, and C. Horgan 

2 It was proved recently that rank-one convexity also implies polyconvexity for isotropic, objective and 
isochoric elastic energies in the two-dimensional case m 
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such as post-buckling and phase transitions etc. Strictly speaking, due to certain necessary 
constitutive constraints such as det F > 0 and objectivity etc, even the domain T a is not 
convex, therefore, it is not appropriate to discuss convexity of the stored energy W(F) in 
general nonlinear elasticity. How to identify global optimal solution has been a fundamental 
challenging problem in nonconvex analysis and computational science. ■ 

Remark 3 (Canonical Duality, Gap Function, and Global Extremality) 

The objectivity is a necessary constraint for any hyper-elastic model. A real-valued function 
W : F a —> M is objective iff there exists a function U (C) such that W(F) = U (F T F) VF E T a - 
By the fact that the right Cauchy-Green tensor C is an objective measure on a convex domain 
£ a = {C E M 3x3 | C = C T , C 0}, it is possible and natural to discuss the convexity of 
17(C). This fact lays a foundation for the canonical duality theory [5], which was developed 
from Gao and Strang’s original work in 1989 [12] for general nonconvex/nonsmooth variational 
problems in finite deformation theory. The key idea of this theory is assuming the existence 
of a geometrically admissible (objective) measure £ = A(F) and a canonical function F(£) 
such that the following canonical transformation holds 

£ = A(F) : T a -> £ a => W(F) = F(A(F)). (10) 

A real-valued function V : £ a — > M is called canonical if the duality relation = VF(£) : £ a —> 
£* is one-to-one and onto [5]- This canonical duality is necessary for modeling natural phe¬ 
nomena. Gao and Strang discovered that the directional derivative A*(F) = 5A(F) is adjoined 
with the equilibrium operator, while its complementary operator A C (F) = A(F)—Aj(F)F leads 
to a so-called complementary gap function, which recovers duality gaps in traditional analysis 
and provides a sufficient condition for identifying both global and local extremal solutions 

mm- ■ 

The canonical duality theory has been applied for solving a large class of nonconvex, 
nonsmooth, discrete problems in multidisciplinary fields of nonlinear analysis, nonconvex 
mechanics, global optimization, and computational sciences, etc. A comprehensive review 
is given recently in m- The main goal of this paper is to show author’s recent analytical 
solutions [7] are valid for general anti-plane shear problems and can be easily generalized for 
solving finite deformation problems governed by generalized neo-Hookean materials. While 
the constitutive constraints in Knowles’ over-determined system p3jj are not necessary, the 
hydrostatic pressure p = p(x i, x 2 ) is independent of X 3 and can’t be considered as a variational 
variable. Some insightful results are obtained on ellipticity condition in nonlinear analysis. 

2 Complete Solutions to Generalized Neo-Hookean Material 

Since the right Cauchy-Green strain C = F 2 F is an objective tensor, its three principal 
invariants 

7 1 (C) = trC, J 2 (C) = ^[(trC) 2 -tr(C) 2 ], J 3 (C) = detC (11) 

are also objective functions of F. Clearly, for isochoric deformations we have / 3 (C) = 1. 
The elastic body is said to be generalized neo-Hookean material if the stored energy depends 
only on I\, i.e. there exists a function V(Ii) such that W(F) = F(/i(C(F))). Since ii = 
tr(F r F) > 0 VF E T a , the domain of V(I\) is a convex (positive) cone £ a = {x E M| x > 0}, 
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it is possible to discuss the convexity of V(Ii). Furthermore, we assume that V(Ii) is a C 2 (£ a ) 


canonical function. Then the canonical transformation (10) for the generalized neo-Hookean 
model is 

£ = A(F)=tr(F T F) : F a -> £ a , W(F) = V(£(¥)). (12) 


For a given external force t(x) on St, we introduce a statically admissible space 
T a = {TGC 1 [B;K 3x3 ]| V-T = 0 in B, N-T = tonSi}. 


(13) 


Thus for any given T G T a , the primal problem (V) for the generalized neo-Hookean material 
can be written in following canonical form 


(V) g : min U( X ) = / P(V*)dH, | V* G * c 


where X c = {x G X a \ A(V%) G £ a } and the integrand P : F a 

P( F) = U(A(F)) — tr(FT). 


is defined by 


(14) 


(15) 


The criticality condition for this canonical variational problem leads to the following 
canonical boundary value problem 

(BVP) g : V • [2£(V*)] = 0 in B, N • [2£(V*)] = t on S t , * = 0 on S x (16) 

which are identical to ( BVP) since VVF(F) = 2F[VU(£)] = 2£F. To solve this fully nonlinear 
boundary value problem is very difficult for direct methods, but easy for the canonical duality 
theory. 

By the canonical assumption of U(£), the duality relation £ = VU(£) : £ a —> £„ is invert¬ 
ible. The complementary energy can be defined uniquely by the Legendre transformation 


^*(0 = {£c-ne)ic = vF(£)}. 


(17) 


Clearly, the function V : £ a —> M is canonical if and only if the following canonical duality 
relations hold on £ a x £* 

C = W(£) £ = VU*(£) o U(£) + U* (£) = ££. (18) 

Using U(£) = ££ — V*(£), the nonconvex function P(F) can be written as the so-called total 
complementary function on X a X £* 

S(F,£) = A(F)£-U*(£)-tr(TF). (19) 

The canonical dual function can be obtained by the canonical dual transformation: 

P d (£) = {~(F,£)| V f ~(F,C) = 0} = —U*(C)-^C-V 2 , r 2 = tr(T T T). (20) 

Thus, the pure complementary energy principle, first proposed in 1998 |3], leads to the fol¬ 
lowing canonical dual variational problem 

(V d ) 9 : sta|n rf (C) = J^P d (C)dB | £ G «S„| , (21) 
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where sta{n rf (£)| C £ S a } stands for finding stationary point of n rf (£) on the canonical dual 
feasible space S a = {£ £ £*\ C ^ 1 " 7-2 £ 

Since the canonical dual variable C is a scalar-valued function, the criticality condition 
<5n d (0 = 0 leads to a so-called canonical dual algebraic equation (see |5j): 

4C 2 W*(C) = t 2 (22) 


Note that W*(C) : £* —> £ a is also one-to-one and onto, this equation has at least one solution 
for any given r 2 = tr(T 7 T) > 0 and C = 0 only if r = 0. Therefore, P d (() is well-defined. 
Due to the nonlinearity, the solution may not be unique [7]. By the pure complementary 
energy principle proposed by Gao in 1999 (see E), we have 

Theorem 1 (Pure Complementary Energy Principle) For any given nontrivial t / 0 
such that T £ T a , pH ) has at least one solution Ck 0, the deformation vector defined by 

Xfc( x ) = \ f C' T ■ fix (23) 

z */xo 

along any path from xq £ S x to x £ B is a critical point ofH(x) and n(x fc ) = n d (£fc). 


This principle shows that by the canonical dual transformation, the nonlinear partial 
differential equation in ( BVP ) for generalized neo-Hookean model can be converted to an 
algebraic equation ( 22 ), which can be solved to obtain a complete set of solutions (see miB!). 

Since S a is nonconvex, in order to identify global and local optimal solutions, we need the 
following convex subsets 


= {C € Sal C > 0}, S+ = {( £ Sal C > 0}, S-= {C € Sale < 0}. (24) 


Then by the canonical duality-triality theory developed in [5] we have the following theorem. 


Theorem 2 Suppose that V : £ a 


is convex and for a given T £ T a such that {£*,} is a 


solution set to (22), F/. = l T, and Xk defined by (23), we have the following statements. 


1. If Ck £ ‘5 Q t_ , then V 2 W(Fk) P 0 and Xk a global minimal solution to (V) g 

2. If fk £ Sf and V 2 lT(Ffc) 0 , then Xk a local minimal solution to (V) g . 

3. If Ck £ and V 2 4T(Ffc) -< 0 , then Xk a local maximal solution to (V) g . 


If {Cfc} C S+, then {xj is a convex set. 
{Ck}cS+. 


The problem (V) g has a unique solution if 


Proof. By using chain rule for W(F) = F(£(F)) we have VTT(F) = 2F[VF(£)] = 2£F, and 

V 2 1T(F) = 2(1 <g> I + 4/i(£)F <g> F, (25) 

where I is an identity tensor in M 3x3 , h(Q = V 2 F(£) > 0 due to the convexity of V on £ a . 
Therefore, V 2 fT(F fc ) y 0 if Ck € S+. 

To prove Xk is a global minimizer of ( (P ), we follow Gao and Strang’s work in 1989 | 12 j . 
By the convexity of V (() on its convex domain £ a , we have 

V(C) - v(Ck) >(C~ Ck)Ck ve, Ck € £ a , Ck = VP( 4 ). (26) 
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For any given variation 5x , we let x = Xk + $X- Then we have 

A(Vx) = tr[(V X ) T (V X )] = A(V Xfc ) + A t (Vx k )(V8x) ~ A c (V<5x), 


( 27 ) 


where A t {F)5F = 2tr[F T (5F)] and A c (5*) = -A(<5*). Clearly, A(F) = A t (F)F + A C (F). 
Then combining the inequality (26) and (27), we have 

n(x)-n(x fe )> f 2C fc tr[(V Xfc ) T (V5x)]d B- f 5 X -tdS+ [ atr[(Vx) T (V*)]d B 

JB JSt JB 


= / [2C k (V Xk ) ~ T] : (V8 X )dB + G ap (5 X , Cfc) V*, G * c 

JB 


(28) 


for any given T G 7)), where 

G ap {x,Q= [ -Ac(V X )CdB= [ Ctr[(V X ) T (V X )]d B (29) 

JB JB 

j. If Xk is a critical point of n(%), then 

[ [2 (V Xk )Ck ~ T] : (V5 X )dB = 0 Vd* G T c , VT G T a 
JB 


is the complementary gap function introduced in 
we have 


Thus, we have II(x) - n(x fc ) > G ap (8x,Ck) > 0 VAy G X c if ( k G S+. This shows that Xk is 
a global minimize!' of {V) g . 

™ i?. u, Io71\ m. T7, — 

2 Sfc 

(30) 


To prove the local extremality, we replace F& in (25) by Ffc = X T such that 


G(Cfc) = v VF(Ffc) = 241 0 I + Q Kik) T 0 T, 

where = W*(Cfc)- Clearly, for a given TgT s such that Cfc G «S“, the Hessian V 2 VF(Ffc) 
could be either positive or negative definite. The total potential n(x fe ) is locally convex if 
the Legendre condition V 2 TT(Vx fc ) >z 0 holds, locally concave if A7 2 W (A7 x k ) ^ 0- Since Xk 
is a global minimizer when Cfc G S(f, therefore, for Cfc G the stationary solution is 
a local minimizer if V 2 IT(Vx fc ) >; 0 and, by the triality theory[5l HI] . Xk is a biggest local 
maximizer if V 2 7F(Vx fc ) -< 0 . 

If {Cfc} C <S^“, then all the solutions {x k } are global minimizers and form a convex set. 
Since P d ( C) is strictly concave on the open convex set <S+, the condition {Cfc} C Sff implies 
the unique solution of (22). In this case, both problems (V) and (BVP) have at most one 
solution. □ 


Theorem 3 (Triality Theory) For any given T G 7 ' a , let Cfc be a critical point of (V° 
the vector Xk be defined by (23), and X Q x S Q C X c x S~ a neighborhood^ of (y fc , Cfc). 

If Cfc G S+, then 


>9> 


X^X 

If Cfc G S~ and G(Cfc) >- 0 , then 


n (Xk) = min n(0 = maxn d (C) = n d (Cfc). 


C6-S + 


n (*fc) = min n(C) = minn d (C) = n d (Cfc). 


If Cfc G 5 a and G(Cfc) -< 0 , f/ien 


n(Xfc) = max n(0 = maxn d (C) = n d (Cfc). 


(31) 

(32) 

(33) 


3 The neighborhood X 0 of Xk * n the canonical duality theory means that Xk i s the only one critical point 
of n(x) on X 0 (see [5]). 
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This theorem shows that the triality theory can be used to identify both global and local 
extremum solutions to the variational problem (V) g and the nonconvex minimum variational 
problem (V) g is canonically equivalent to the following concave maximization problem over 
an open convex set 5+, i.e. 

( 1 P #) 9 : max jn d (C) = J^P d (()dB \ ( E 5+1 , (34) 

which is much easier to solve for obtaining global optimal solution of (V) g . 


3 Generalized Quasiconvexity, G-Ellipticity, and Uniqueness 


The ellipticity is a classical concept originally from linear partial differential systems, where 
the deformation is a scalar-valued function x '■ & —> M and stored energy is a quadratic 
function W( 7) = ^7 T H7 of 7 = Vy E M 3 . The linear operator 

L[ x ] = -V • [H(Vy)] = -[hijXjU 

is called elliptic if H = {hij} is positive definite. In this case, the function P( 7) = W (7)—7 t t 
is convex and its level set is an ellipse for any given r E M 3 . This concept has been extended 
to nonlinear analysis. The fully nonlinear partial differential equation in (BVP ) ([ 5 ]) is called 
elliptic if the stored energy W (F) is rank-one convex. In the case Ik G C 2 , the rank-one 
convexity is equivalent to the Legendre-Hadamard (LH) condition: 

(a <g> a) : V 2 IT(F) : (r? <g> rj) > 0 Va, t] E M 3 , VFeT a C M mxn . ( 35 ) 


The (BVP) is called strong elliptic if the inequality holds strictly. In this case, (BVP) has 
at most one solution. In vector space, the LH condition is equivalent to Legendre condition 
V 2 IF( 7 ) 7 0 V 7 E M n . 

Clearly, the LH condition is only a sufficient condition for local minimizer of the variational 
problem ( V). In order to identify ellipticity, one must to check LH condition for all local 
solutions, which is impossible for general fully nonlinear problems. Also, the traditional 
ellipticity definition depends only on the stored energy W(F) regardless of the linear term in 
P(F). This definition works only for convex systems since the linear term can’t change the 
convexity of P(F). But this is not true for nonconvex systems. To see this, let us consider 
the St. Venant-Kirchhoff material 

W(F) = 1e:H:E, E = ^[(F) T (F) — I], (36) 


where I is a unit tensor in M 3x3 . Clearly, this function is not even rank-one convex. A 
special case of this model is the well-known double-well potential W( 7 ) = |(^| 7| 2 — l) 2 . In 
this case, if we let £ = A( 7 ) = ^| 7| 2 — 1 be an objective measure, we have the canonical 
function V(£) = ^£ 2 . In this case, the canonical dual algebraic equation (22) is a cubic 
equation (see 151) 2 £ 2 (C + 1 ) = t 2 , which has at most three real solutions {£jT(x)} at each 
xgB satisfying £1 > 0 > £2 > £ 3 - It w &s proved in j5] (Theorem 3.4.4, page 133) that for a 
given force t(x), if t 2 (x) > 8/27 Vx £BcK, then (BVP) g has only one solution on B. If 
t 2 (x) < 8/27 Vx E B s C B, then (BVP) g has three solutions (xfc(x)} at each x E B s , i.e. 
n(y) is nonconvex on B s . It was shown by Gao and Ogden that these solutions are nonsmooth 
if t(x) changes its sign on B s m- 
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Analytical solutions for general 3-D finite deformation problem (V) were first proposed 
by Gao in 1998-1999 [31 3]. It is proved recently [ 8 ] that for St Venant-Kirchhoff material, 
the problem (V) could have 24 critical solutions at each material point x G B, but only one 
global minimizer. The solution is unique if the external force is sufficiently bigger. 

The geometrical explanation for ellipticity and Theorem [2] is illustrated by Fig. [lj which 
shows that the nonconvex function P( 7 ) = ^^ |^y| 2 — l) 2 — 7 T r depends sensitively on the 
external force t G R 2 . If |r| is bigger enough, P( 7 ) has only one minimizer and its level set 
is an ellipse (Fig. [l] (b)). Otherwise, P( 7 ) has multiple local minimizers and its level set is 
not an ellipse. For r = 0, it is well-known Mexican-hat in theoretical physics (Fig. 0 (a)). 



Figure 1: Graphs and level sets of -P(x) with r = 0 (left) and t/ 0 (right) 

Fig. a shows that although P( 7 ) has only one global minimizer for certain given r, the 
function is still nonconvex. Such a function is called quasiconvex in the context of global 
optimization. In order to distinguish this type of functions with Morry’s quasiconvexity in 
nonconvex analysis, a generalized definition on a tensor space F a C R mxn could be convenient. 

Definition 1 (G-Quasiconvexity) A function P : T a C R mxn —> R is called G-quasiconvex 
if its domain F a is convex and 

P(0 F + (1 - 0)T) < max{P(F),P(T)} VF, T G F a , V0 G [0,1], (37) 

It is called strictly G-quasiconvex if the inequality holds strictly. 

For a given function P : F a —> R, its level set and sub-level set of height a G R are defined, 
respectively, as the following 

C a (P) = {F G Fa P( F) = a}, &(P) = {FeF a \ P(F) < a}, a G R. (38) 

Moreover, we may need a generalize ellipticity definition for nonconvex systems. 

Definition 2 (G-Ellipticity) For a given function P : F a —> M and a G M, its level set 
C a (P) is said to be a G-ellipse if it is a closed, simply connected set. For a given t such that 
T G T a , the ( BVP ) is said to be G-elliptic if the total potential function -P(F) is G-quasiconvex 
on F a . (BVP) is strongly elliptic if P( F) is strictly G-quasiconvex. 

Clearly, we have the following statements: 

P(F) is G-quasiconvex & a (P) is convex C a (P) is G-ellipse Va G R . (39) 

P(F) convex => rank-one convex => G-quasiconvex => (BVP) is G-elliptic. (40) 
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This statement shows a fact in nonconvex systems, i.e. the number of solutions to a 
nonlinear equation depends not only on the stored energy, but also (mainly) on the external 
force field. The nonlinear partial differential equation in ( BVP ) is elliptic only if it is G- 
elliptic. (BVP) has at most one solution if the integrand P( F) in the total potential n(%) is 
strictly G-quasiconvex on T a . 

In global optimization, the most simple quadratic integer programming problem 

(V)i : min |ll(x) = -x 1 Qx — x T t | x £ {0, l} n C M n 

could have up to 2 n local minimizers due to the indefinite matrix Q and the integer constraint. 
Such a nonconvex discrete optimization problem is considered as NP-hard in computer science. 
However, by using canonical transformation £ = {xi(xi— 1)}, the canonical dual of this discrete 
problem is a concave maximization over a convex set in continuous space. It was proved in | 6 ] 
that as long as the source term t = {tj} is bigger enough, 5+ / 0 and (V)i is not NP-hard. 
The decision variable is simply {xi} = {0 if £j < 0, 1 if ti > 0} (Theorem 8 , [ 6 ]). 


4 Anti-plane Shear Deformation Problems 


Now let us consider a special case that the homogeneous elastic body B C M 3 is a cylinder 
with generators parallel to the e 3 axis and with cross section a sufficiently nice region O C M 2 
in the ei x e 2 plane. The so-called anti-plane shear deformation is defined by (see Knowles 
1976, [I!]) 

X(x) = |a - 2£i, A _ 5x 2 , Ax 3 + u(x\, x 2 )j : Fl —> M 3 , (41) 

where (aq, x 2 , x%) are cylindrical coordinates in the reference configuration B relative to a 
cylindrical basis {e*}, i = 1,2,3, the parameter A is a positive constant, and u : Fl —> M is 
the amount of shear (locally a simple shear) in the planes normal to e 3 . On T x C <9P, the 
homogenous boundary condition is given u(x a ) = 0 Vx a £ F v , a = 1,2. On the remaining 
boundary T t = <90 n T x , the cylinder is subjected to the shear force 


t(x) = f(x)e 3 Vx e F t , 


where t : —>• M is a prescribed function. According to Knowles, the deformation gg may 
be thought of as one in which the body first undergoes an axial elongation (or contraction) 
of stretch ratio A (regarded as given), and is then subjected to an anti-plane shear with 
out-of-plane displacement u. For this anti-plane shear deformation we have 


( \-\ 

o 

o 

( A- 1 + nl 

U,lU }2 

A u t i \ 

0 

A“2 o 

, c = F t F = U 1U , 2 

A -1 + u 2 2 

A U t 2 

V «,i 

U ,2 A ) 

V Aur 

A U y 2 

A 2 J 


where u ;Q represents du/dx a for a = 1,2. By the notation |Vu | 2 = u 2 x + u 2 2 , we have 
/ 1 (C) = A 1 + |Vu| 2 , / 2 (C) = A 2 + A- 1 |Vu| 2 , J 3 (C) = 1 , 
where Ai = A 2 + 2A _1 , A 2 = A -2 + 2A. Particularly, 

/i(C) = / 2 (C) = 3 + |Vu | 2 if A = 1. 


(42) 


(43) 

(44) 
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Lemma 1 For any given A > 0, the homogenous hyper-elasticity for general anti-plane shear 
deformation must be governed by a generalized neo-Hookean model, i.e. W(F) = V(I\). 
Particularly, the Mooney-Rivlin model is identical to the neo-Hookean model subjected to a 
constant, i.e. W(F) = A{I\ — 3) + B with B = 0 if A = 1. 

The proof is elementary, i.e. by the fact that I 2 = A _1 /i + a, a = X 2 — A - 1 Ai, we have 
W( F) = W(h,I 2 ) = W(h, A" 1 /! + a) = V (ii), VA > 0. 

Moreover, let A = c\ + C 2 A _1 , B = C 2 ( 3 A -1 — 3 + A 2 — A 1 A -1 ), we have the Mooney-Rivlin 
model W{ F) = ci(7i - 3) + c 2 (/ 2 - 3) = A{h - 3) + B. 


The fact detF = 1 shows that the anti-plane shear state (41) is an isochoric deformation. 
Therefore, the kinetically admissible displacement space X c can be simply replaced by a 
convex set 

U c = {u(xi,X 2 ) £ M)| u(x a ) = 0 \/x a £ T x }. (45) 

Thus, in terms of 7 = Vu, f = A( 7 ) = I 71 2 + Ai, and VU(F( 7 )) = V(A(-f)), for any given 

T6Ta = {TGC 1 [H;K 2 ]| V-r = 0 in J7, n r = t on R} 

Problem ( V) g for the general anti-plane shear deformation problem has the following form 


( V)s : min |n(u) = f P(Vu)dtt \ u E W c | , P{ 7 ) = V(A( 7 )) - 7 T r (46) 
Under certain regularity conditions, the associated mixed boundary value problem is 


(BVP) S : 


( V • (2CVii) = 0 in U, 

\ n • ( 2 CVrt) = t on R, u = 0 on T x 


(47) 


where n = {n a } £ M 2 is a unit vector norm to <9U, and £ = VF(£), | = Ai + |Vit|. 

If T x = dfl, then ( BVP) S is a Dirichlet boundary value problem, which has only trivial 
solution due to zero input. For Neumann boundary value problem Tj = dQ, the external 
force field must be such that 


/ t(x)dT = 0 

Jr t 


for overall force equilibrium. In this case, if x is a solution to ( BVP ) S , then x = X + c is 
also a solution for any vector c £ M 3 since the cylinder is not fixed. 

By the fact that the only unknown u is a scalar-valued function, anti-plane shear defor¬ 
mations are one of the simplest classes of deformations that solids can undergo m- Indeed, 
if V(£) is a canonical function and for any given r £ F a such that r = |r|, the canonical dual 
problem has a very simple form 


(V d ) s : sta|n d (C ) = jT 


AiC-^O-^CV 


dfl | C £ Sa 


(48) 


Since A (it) = \Vu\ 2 + Ai, the canonical dual equation ( 22 ) for this problem is 


4C 2 [VU*(C) - Ai] = r 2 . 


(49) 
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Theorem 4 For any given pre-stretch A > 0 and non-trivial shear force f(x) 7 I 0 onTt such 
that T a / 0, the canonical dual problem ( V d ) s has at least one non-trivial solution and 


1 f x 

M*) = 9 / Cfe 1 

z -'xo 


r • dx 


(50) 


a/ong any path from xq £ to x e fi b a critical point ofU(u) and Il(ufc) = LI rf (£*.). 


Proof. By the fact that the nontrivial shear force t leads to t > 0. The equation (49) has 
at least one nontrivial solution 7 I 0. By the pure complementary energy principle, the 
equation (50) gives a nontrivial solution to the anti-plane shear deformation. □ 


Note that F is an affine function of 7 = Vu G M 2 , it is also mathematically equivalent to 
assume the existence of a real-valued function W : M 2 — > M such that 

W(F(u)) = V(ii(Vu)) = W{j(u)) Vu € U c (51) 

holds for general anti-plane shear deformation problems without any additional constitutive 
constraints. In this case, by choosing £ = A(Vtt) = |Vu| 2 and the canonical transformation 
^(£( 7 )) = W"(7), equivalent results for complete set of solutions have been obtained for both 
convex and nonconvex anti-plane shear deformation problems |7j. 

Clearly, the anti-plane shear problem ( BVP) S is linear only if V{I\) = A(I\ — 3) for a 
given constant A > 0, i.e. the neo-Hookean model. For nonlinear elasticity, the problem (V) s 
could have multiple critical solutions {iifc( x )} at each x £ fl s C fi. As long as 7 ^ 0 , the 
boundary value problem ( BVP) S should have infinitely many solutions (see [10]). Therefore, 
it is impossible to use Legendre condition to identify global minimal solution. Theorem [2] 
shows that the Legendre condition is only necessary but not sufficient condition for global 
optimality. The sufficient condition is simply 

Ck £ <5^ 44 Gap = [ Cfc|Vu| 2 dfl > 0 Vu G U c , Cfc £ , (52) 

in 

which was first proposed in 1992 0 . If all solutions {Cfc} £ <S+, Problem (T) s is G- 
quasiconvex, which has unique solution if {Cfc} £ ^q 1 ". Application of Theorem [4] has been 
illustrated for both convex and nonconvex problems given recently in [7|. 


5 Remarks on Knowles’ over determined problem 


Now let us revisit Knowles’ work in 1976 nn- Instead of the minimal potential variational 
problem ( V ), Knowles started from the strong form of (V), i.e. diver = 0 in the boundary 
value problem ( BVP) p given in ([9]) with general constitutive law for incompressible materials 


_ dw(h,i 2 ) 

a dF 


— pF~ T . 


(53) 


For the same anti-plane shear deformation problem (41), he ended up with three equilibrium 
equations (i.e. equations (2.19) and (2.20) in 11 11 


q,a + {2W2U j0l U t p) p - p,3U,a = 0, 

4 There is a mistake in M, i-e. W\ in Knowles’ equation (2.19) should be W 2 


(54) 
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(55) 


[2 {Wt + ^ - A-^s = 0, 

where W a = dW/dI Q = ( a , a, /3 = 1, 2 and q = Xp - 2m - 2(A 2 + A " 1 + |Vu| 2 )IK 2 . 


The hrst two equations in (54) are corresponding to the general equilibrium equation 


= 0 in ei and e 2 directions; while the third one (55) is in e 3 direction. Knowles indicated 


(Equation (2.22) in [Tl] ) that the hydrostatic pressure p is linear in X 3 , i.e. 

p = cx 3 +p(x i,x 2 ) 


(56) 


where c is a constant. Saccomandi emphasized recently that p is the Lagrange multiple 


associated with the incompressibility constraint, which must be in the form of (56) and c / 0 
for general incompressible material 


Clearly, for a given strain energy IK(F) = W(I 1 , 12 ), the governing equations obtained by 
Knowles constitute an over-determined system in general, i.e. two unknowns (u,p) but three 
equations. In order to solve this over determined problem, Knowles believed that the stored 
energy W(R , I 2 ) should have some restrictions and he proved the following theorem. 

Theorem (Knowles, 1976 |14]) If the stored energy W(R,l 2 ) is such that the ellipticity 
condition (i.e. the equation (3.5) in m 

d.[2R{W\ + A -1 ^)] 


dR 


>0 VR > 0, A > 0 


(57) 


holds, then the associated incompressible elastic material admits nontrivial states of anti-plane 
shear for a given pre-stretch A if and only if W(R, I 2 ) also satisfies the following constitutive 
constraint (i.e. equation (3.22) in m 

bm + ( 6 A _1 - 1)W 2 = 0 , (58) 

for some constant b, for all values of R, I 2 such that I\ = A 1 + -R 2 , I 2 = A 2 + A _ 1 i? 2 , R = |^y|. 

First, by Lemma[l]we know that W(F) = W(I\, I 2 ) = V(I\) hold for any given anti-plane 
shear deformation. There is no need to have both R, R as variables. Therefore, the following 


trivial result shows immediately that Knowles’ condition (58) is not a constitutive constraint. 


Lemma 2 For any given stored energy W(F) = W(R,I 2 ) such that R = Ai + |Vu| 2 , I 2 = 
X 2 + A _ 1 |Vu| 2 , Knowles’ constitutive condition (58) is automatically satisfied for b = ^A. 


The proof of this statement is elementary: by chain rule and Ii = XR + Ai — AA 2 , we have 
immediately 


- dWdR 
W 2 = wy-ww = XWi 
oil 0 I 2 


bWi + ( 6 A _1 - 1 )W 2 = (2b - X)W\ = 0 Vb = f-X. 


1 


To check if the Lagrange multiplier p = p(xi, X 2 , X 3 ) must be in Knowles’ formula (56), 
we use mathematical theory of Lagrange duality. For any given real-valued function </>(x) E 
L q (B ), the Lagrange multiplier p for the equality constraint </>(x) = 0 must be in the dual 
space L q (£>) such that 1/q + 1/q' = 1. Since F depends only on (x\,X 2 ) 6 12, the constraint 
^>(x) = detF(rt) — 1 is defined on 12 C B, its Lagrange multiplier p(x) must be defined on 12. 


Indeed, by simple calculation for the form (56) 


(t)(xi,X2)p(xi,X2,X3)&B = / f>(x 1 ,X2)[ p(xi,X2,x 3 )dx 3 }dQ 

: Jq J 
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one can easily find that the Lagrange multiplier is independent of X 3 . Thus, we must have 
c = 0 and p = p(x 1 , 2 : 2 ) for any anti-plane shear deformation. For this reason and IV 2 = AVFi, 
c = W(£) = Wi, the equation (55) (i.e. (2.20) in (T3J) is identical to the equation in ( BVP ) S : 


[2(Wi + \~ 1 W 2 )ufi\ 3 = 0 V • [CVu] = 0 in U. 


(59) 


Now we need to check the other two equilibrium equations in Knowles’ over-determined 
system. Instead of the local analysis, we use the well-known virtual work principle 


[ tr(a ■ 8F(x))dB = I t • 5x, € X ( 

JB JS t 


(60) 


which holds for any given deformation problem regardless of constitutive laws. For smooth de¬ 
formation x and sufficiently regular B and dB, we have the following strong complementarity 
conditions 

(<5x) • (diver) = 0 in B , (d'x) ■ cr ■ N = (8x) ■ t on S t (61) 


The fact that the anti-plane shear deformation (41) has no displacements {rq} in ei and e 2 
directions, i.e. 5\a = 0 Va = 1,2 a.e. in 11, the vector diver is not necessarily to be zero in 
these directions. This shows that the additional two equilibrium equations (54), i.e. (2.19) 
in the paper m , can’t be obtained from the virtual work principle. By the fact that the 
boundary value problem ( BVP) S is well-determined by the equation (59), these two extra 
equations are useless for the problem considered. 

To understand the “function” of the hydrostatic pressure p(x) in Knowles’ over-determined 
problems for either compressible or incompressible materials, we use the KKT complemen¬ 
tarity condition in ([ 8 ]), i.e. p(detF — 1) = 0. As we know that the anti-plane shear state is a 
volume preserving deformation, the equality detF(ri) = 1 is trivially satisfied all most every 
where in Q for any materials. Thus, we must have p(x) ^ 0 a.e. in 14, i.e. the only function 


of this arbitrary non zero parameter is to balance the extra two equations (54), which can’t 


be obtained by the virtual work principle. This shows that the governing equations obtained 
by the minimum total potential principle are always compatible. 

Finally, let us exam the ellipticity condition in Knowles’s theorem. On page 407 of m, 


Knowles indicated: the condition (57) “guarantees that (59) is elliptic at every solution u and 
at every point in 0”. The following theorem is important in nonlinear analysis. 


Theorem 5 The ellipticity condition (57) is neither necessary nor sufficient for the nonlinear 
PDE [59 1 ) to admit nontrivial states of anti-plane shear. ( BVP) S has at least one solution 
only if t(x) / 0 onT t such that T a / 0. 

For any given convex function W(I\, I 2 ) and the external force t(x) /0 onT t , the equation 


(59) is strongly G-elliptic if 


Ci > 0, Ci = Wi(/i,/ 2 ) 


(62) 


for every solution Ci of 


Proof. Let £ = By using chain rule for W^) = 1F(C(7)) 

W(7) = V 7 TK(C(7)) = 2 7 (TKi + A” 1 !^), 


thus, Knowles’ ellipticity condition ([57]) is actually a special case of the strong Legendre 
condition V 2 1K( 7 ) >- 0, which can only guarantee the convexity of VF(F) = 1F(7), i.e. under 
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this condition, the ( BVP) S has at most one solution. Clarly, ( BVP) S has a trivial solution 


if f(x) = 0 on iy Therefore, Knowles’ ellipticity condition (57) is not sufficient to admit a 
nontrivial solution. 

By the canonical duality theory we know that for nonconvex stored energy W{ F) = W ( 7 ), 
the (BVP) S has multiple nontrivial solutions if t(x) ^ 0 on K such that T a 7 ^ 0 - Therefore, 


Knowles’ ellipticity condition (57) is also not necessary to admit a nontrivial solution. 
By simple calculation for ( |57[ ), we have 

2 (Wi + \- 1 W 2 ) + 4R 2 (W n + 2X~ 1 Wu + \~ 2 W 22 ) > 0 , 


(63) 


which is a strong case for (25), where W a p = d 2 W/dl a dlp. If the canonical function W(I\, I 2 ) 
is convex in £ = {I \, / 2 }, we have 

Wu + 2 A- 1 IKi 2 + A- 2 IK 22 > 0 V{/i,/ 2 } e« 2 , A > 0 . (64) 

By the facts that ( 2 = W 2 = XW± = A(j and £1 = VK (I\) = W\ , we know that the condition 


(63) holds as long as 

2(W\ + A _ 1 IK 2 ) = 4Ci > 0. 

Thus, by Theorem [ 2 ] we know that the function P( 7 ) is strictly G-quasiconvex and (59) is 
strongly G-elliptic. In this case, ( BVP) S has at most one solution. □ 

Combining Theorems m and Lemma [2] we know that Knowles’ constitutive constraints 
(57) and (58) are neither necessary nor sufficient for the existence of nontrivial states of anti¬ 


plane shear. Actually, this ellipticity condition even disallows many possible nontrivial local 
solutions in nonconvex problems. Indeed, it was shown in mu that for any given nonconvex 
stored energy IK(F( 7 )) = W(hij) ,/ 2 ( 7 )) = IK( 7 ) and nontrivial external force f(x) / 0 , 
the minimum potential variational problem ( V) s has at least one solution {uk} in Banach 
space U c , which can be obtained analytically by the canonical duality theory. If f(x) is very 
small, the solution may not unique, the one such that Ci u k) £ <5^ is a global minimal solution. 
Both global and local minimum solutions could be nonsmooth if r(x) changes its sign in fh 
While Knowles’ over-determined system admits only a unique smooth solution in C 2 due to 
the additional ellipticity restriction on W. Therefore, Knowles’ over-determined system is a 
very special case of the variational problem ( V) s . 


6 Conclusions 

In summary, the following conclusions can be obtained. 

1. The pure complementary energy principle and canonical duality-triality the¬ 
ory developed in [5j are useful for solving general nonlinear boundary value 
problems in nonlinear elasticity. 

2. The ellipticity condition for fully nonlinear boundary value problems in finite 
deformation theory depends not only on the stored energy function, but also 
on the external force field. 

3. The triality theory provides a sufficient condition to identify both global and 
local extremum solutions for nonconvex problems. 

4. General anti-plane shear deformation problems must be governed by the gen¬ 
eralized neo-Hookean model. 
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5. Unless the KKT theory is wrong, the incompressibility is not a variational 
constraint for any anti-plane shear deformation problem, the pseudo-Lagrange 
multiplier p depends only on (x\,X 2 ), which is not a variable for the problem. 

6 . Unless the virtual work principle is wrong, there is only one equilibrium equa¬ 
tion for general anti-plane shear deformation problems. The two extra equa¬ 
tions in Knowles’ over-determined system are not required. 

7. Unless the minimum potential variational principle is wrong, the constitutive 
conditions required by Knowles’ Theorems in mm are neither necessary 
nor sufficient for general homogeneous materials to admit nontrivial states of 
anti-plane shear. 

The first three conclusions are naturally included in the canonical duality-triality theory 
developed by the author and his co-workers during the last 25 years [5] . Extensive applications 
have been given in multidisciplinary fields of biology, chaotic dynamics, computational me¬ 
chanics, information theory, phase transitions, post-buckling, operations research, industrial 
and systems engineering, etc. (see recent review article EH)- 

The last four conclusions are obtained recently when the author got involved in the dis¬ 
cussions with colleagues on anti-plane shear deformation problems. As highly cited papers 
mm, Knowles’ over-determined system has been extensively applied to many anti-plane 
shear deformation problems in literature, see recent papers PU20U2I1E2]. This is the moti¬ 
vation for this paper. 
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